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ABSTRACT 

Dynamic equations which include the effects of 
unsteady aerodynamic forces and a flexible body 
structure have been developed for a free-flying 
high performance fighter aircraft. The linear and 
angular deformations are assumed to be small in 
the body reference frame, allowing the equations 
to be linearized in the deformation variables* 
Equations for total body dynamics and flexible 
body dynamics are formulated using the hybrid 
coordinate method and integrated in a state space 
format. A detailed finite element model of a 
generic high-performance fighter aircraft is used 
to generate the mass and stiffness matrices. 
Unsteady aerodynamics are represented by a 
rational function approximation of the doublet 
lattice matrices. The equations simplify for the 
case of constant angular rate of the body 
reference frame, allowing the effect of roll rate 
to be studied by computing the eigenvalues of the 
system. It is found that the rigid body modes of 
the aircraft are greatly affected by introducing a 
constant roll rate, while the effect on the 
flexible modes is minimal for this configuration. 

1. INTRODUCTION 

Future fighter aircraft must be able to meet 
stringent maneuverability and performance 
requirements. This will result in aircraft 
designs in which the interaction of flexibility, 
aerodynamics, and overall body motion during a 
maneuver are of prime importance. The need for 
superagility and the use of advanced lightweight 
materials will make it very important to consider 
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flexibility effects in the analysis of the 
aircraft undergoing maneuvers at high rates. 

Flexible body dynamics have been investigated 
in many other writings, including references [1] - 
[ 4 J . In this paper, dynamic equations will be 
derived in a manner similar to that in reference 
[5], which contains a more thorough development of 
the equations. In addition, aerodynamic forces 
will be explicitly included in the equations. 
These equations will then be applied to a 
realistic model of a modern fighter aircraft. 

The aircraft is assumed to be a collection of 
elastically interconnected, discrete rigid 
subbodies which are subjected to external forces 
and torques, including unsteady aerodynamic 
forces. It is assumed that the deformations of 
the subbodies with respect to the body reference 
frame are small so that the high order terms in 
the deformation variables and their rates can be 
neglected. The rotational effects of motors, 
fans, and turbines are not included in this 
representation. 


BODY REFERENCE FRAME 



Figure 1. Reference Frames and a Subbody 


145 



Table 1. Vectors and Dyads definitions 


Vectors and Dyads 

Definition 

(i)T 

( b )I 

{ i J T = {b) T 9 
« - (b} T (o 

Tb s ) T 

d = {i)J d 

V - (i}T V 

Is - (U T *s 

c . {b) T c 
£s - (bl* r s 
Hs - ( b ) T - u s 

is - (M T p s 

is - (bl! P s 
I - {b)T I (b } 

Is " f b s ) I s {b s ) 

basis matrix for inertial reference frame (IRF) 

basis matrix for body reference frame (BRF) 

direction cosine matrix relationship between IRF and BRF 

angular velocity of body reference frame 

basis matrix for sth elemental body reference frame(ERF) 

Position of current center of mass (CM) in IRF 

Velocity of current CM in IRF 

Position of the sth element in IRF 

position of CM in BRF 

position of undeformed sth element from undeformed CM 
position of sth element from the undeformed position 
position of differential mass in sth elemental body 
angular deformation of sth element 
Inertia dyadic of the aircraft with respect to CM 
Inertia dyadic of the sth elemental body about its CM 


Table 2. Vector Identities and Matrix Operation 
Equivalents 


Vector representation: 


(i) 1 


Z « Ci } 1 


x 1 
x 2 

Vi 

Y2 

^3- 


- U ) 


= {i } 1 


Cross product representation: 

X X l = {i) T X Y = - { i ) T Y X 
where the ~ operator is defined below: 


<xr = x = 


0 -X 3 x 2 

x 3 0 ~ X 1 

L-x 2 xj 0 


Operations with dyad 

I - (i} T I UJ 
I * x « { i } T I X 
l x x • i = ( i } T X X Y 


A highly detailed description of the hybrid 
coordinate method, which is used here to develop 
the dynamic equations, can be found in references 
[5] and [6], Only some highlights of the 
development of these equations will be presented 
in this paper. The development of the equations 
closely follows that of reference [7], with 
aerodynamic forces added. The equations are 
implemented as a computer program, FLXAIR, 


Figure 1 shows a schematic diagram of the 


various reference frames associated with each 
subbody. Definitions of the vectors and dyads 
used in this figure and in the derivations are 
given in Table 1. Table 2 shows the equivalency 
between various operations in a vector/dyad format 
and those in a matrix format. The matrix format 
is used for implementing the computer solution to 
our problem. 

Section 2 deals with the net force and torque 
applied to the total body. The net forces and 
torques on the subbodies are described in section 
3. Derivations are kept brief, with only main 
steps provided. The rational function 

approximation used for describing the unsteady 
aerodynamic forces Is given in section 4. In 
section 5 the equations are integrated in a state 
space format, with aerodynamic forces specifically 
separated from other external forces. These 
nonlinear and time dependent equations can be used 
for simulation. Vhen the angular velocity of the 
body reference frame is constant, the equations 
become time invariant. It is then possible to 
study the effects of angular velocity on vehicle 
structural dynamics by performing an eigenvalue 
analysis. Vhen unsteady aerodynamic loading is 
included in this formulation, this is seen to be a 
flutter analysis under maneuver. 

A large-order finite element model which is a 
realistic representation of an advanced fighter 
was used to demonstrate the stability effects of 
high roll rates. Section 6 describes the NASTRAN 
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model, structural and aerodynamic, which vas used 
in the analysis. Results are shown for various 
roll rates and for variations of overall stiffness 
of the aircraft. The analyses show little effect 
on the flexible modes of the system due to roll 
maneuvers. Considerable effect was, however, 
observed for the rigid body modes. 

2. TOTAL BODT DYNAMICS 

The equations are derived from Newton-Euler 
equations. The equations for the net force F, and 
the net torque P , can be represented as follows. 

F . ^"Tot d) (1) 

P = (H) (2) 

dt 

M Tot is the total mass, and H is the angular 
momentum referred to the CM of the aircraft. 
Presuperscr ipt i refers to the fact that the 
differentiation must be with respect to the 

inertial reference frame. 

The further development neglects the effects 
of rotating bodies such as engine compressors, 
fans, rotors, etc.. It is assumed that the 

deformation of flexible bodies is small in the 
body reference frame. This assumption is used to 
neglect the high order terms in the deformation 
variables u s (linear deformation of sth element) 
and 3 s (angular deformation of the sth element) and 
their derivatives. 

Equation (1) can be written in the body 
reference frame and for ease of computer 
implementation in matrix form as follows: 

F - M Tot 0V (3) 

The development of equation (2) to a computer 
implementable step is lengthy. Only few key steps 
are given. 

The angular momentum H is defined as 

H - /(£ 4 r s +u s + p s ) x *d (c+r 5 +u s +p s ) dm (4) 

3T 

The development makes use of the mass-center 
definition 

/(c+r^Ug+pj) dm - 0 (5) 


and the following identity: 

f(£s + Hs*P s> x ( “ * (£s*Hs + £s) ) dm - I •« (6) 

With the use of equations (4) to (6), equation (2) 
can be written as 

P = I • u> + w x X ‘ w + X*w + c x c 

+ l 6 J(r s +u s +p s ) x (£ s 4U s +p s ) dm (7) 

d t 

The assumptions of discrete lumped masses and 
small deformations with respect to the body 
reference frame are now used to convert the 
integration operation into the following summation 
operation : 

f(£s + Ms + £s) x (£s*Ms + £s) m 

l ls x m s Hs + E Js ' < 8 > 

where m s is defined as: 

J dm = m s (9) 

s 

and center of mass definition of sth lumped mass 
is given by 

J p s dm = 0 (10) 

s 

Finally, equation (7) can be written in a 
computer implementable form as: 

P = Iw +Iw +wlu> 

+ I(((wr s )“ +r s u>)u s m s +r s u s m s 4l s S s ) 

4 «n s e s (11) 

The total inertia is assumed to be linear in 
the deformation variables. 

I - I*4Em s (2r s T u s E -r s u| -u s r s T ) 

4 i(a s i s -i s B s ) < 12 ) 

where I* is the inertia of the undeformed 

airplane. Therefore 

i . £m s (2r s T i s E -r s GJ -u s r s T ) 

. r(e s r s -i s e s ) (13) 

Equations (12) and (13) can be substituted into 

equation (11) to further simplify the equation. 
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3. ELEMENTAL BOOT DYNAMICS 


4. RATIONAL FUNCTION APPROXIMATION 
OP UNSTEADY AERODYNAMICS 


The net forces and torques on the sth 

elemental body are as follows: 

is - m s * £ ♦ Is + Hs> < 14 > 

at z 

2s - ^ H s (15) 

dt 

H s is the inertial angular momentum of sth element 
referred to its mass center. Note that this 
equation is applicable to all n subbodies. 

It is assumed that the body reference frame 
and the elemental body reference frames are 
initially colinear. This assumption, though not 
necessary, is used here to simplify the equations. 

Equation (14) can be written in the following 
computer implemen table form: 

f s = m s { 0V + c+ w(r s +u s +c) +2w{c+u s ) 

+ u s +to£>(r s + u s + c) ) (16) 

noting that 

V = d_(d) (17) 

dt 

H s is defined as 

H s - Is • « (18) 

Invoking the assumption of small deformation, the 
rotation is represented by 

e s - (b) T @ s (19) 

Note that this equation is strictly true if the 
rotations are infinitesimally small. The 
relationship betveen the body reference frame and 
the elemental reference frame can now be 
approximated as 

(b s ] T - (b) T (E+£ s ) (20) 

where E is a 3x3 unit matrix. 

Using equations (18) to (20), equation (15) 
can be written in the following computer 
iroplemen table form: 

Ps " 1 3 ( <*>* P 5 ) ♦(I s k+«l s -<l s u> r )0 5 +S>I s W 

♦ (I s S) -(I s w)“ -«{I s m)" ♦5>I s «)0 s (21) 


The formulation of the unsteady aerodynamics 
is based on the relation 

( 7 ) - -K INIDl (Ap) (22) 

P» 

where (Ap) represents the pressures at aerodynamic 
force nodes, (v) contains the velocities normal to 
the lifting surface induced by {Ap}, and [NIDJ is 
the induced normal dovnvash influence matrix. The 
induced velocities are defined as dovnvash 
collocation points which are located at the 3/4 
chord of each aerodynamic box for the doublet 
lattice method. 

Dovnvash collocation points are those points 
on a lifting surface at which the induced velocity 
normalized by the free stream velocity is equal to 
the local angle of attack (a), i.e., 

(a) - ( 7 ) (23) 

The pressures are then given by 

(Ap) - j pV 2 [NIDI " 1 (a) (24) 

or 

(bp) - \ pV 2 |AIC] (a) (25) 

where [AIC] - I NXD J — 1 . 

In the following derivation, Equation (24) is 
used as a starting point. 

From (Ap), by an integration or "lumping" 
process represented by (ZF), the aerodynamic 
forces are obtained: 

(Z aero ) - [ZF] (Ap) (26) 

The local angle of attack, taken relative to 
the free stream velocity, V, is given by 

(a) - (a^.) ♦ (a-) (27) 

The contribution a T is the instantaneous 
slope of the lifting surface, relative to V, in a 
plane through V perpendicular to the lifting 
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surface: 


and 


U T ) - ID 0 1 [i} (28a) 

where [D 0 J is a differentiating matrix. 

The contribution a z results from the rate of 
translation in a direction perpendicular to the 
lifting surface: 

l*i» - 0} (20b) 

where l D z ] is an interpolating matrix. 

Substituting Equation (28) into Equation (27) 
and replacing z by sz yields: 

(a) = [l D q] ♦ | |D,]] (z) (29) 

Combining Equations (24), (26) and (29) leads 
to : 

(Zaero) - \ [NID]' 1 

* [l D 0 ] * | (D,]]{Z) (30) 

For constant amplitude oscillation s = i(0 - 
i(Vk/c). The induced velocity matrix is a 
function of ik. It follows that Equation (30) can 
be written as: 


|A z (lk)l - |Zp] INID(lk)!- 1 |D Z ] 


{Z aero ) - | o'* 2 lAfik.)] (z) 


(31) 


where A(ik) is given by: 

[A(lk)] - [ZP] [NID(ik)!- 1 [(D 0 ) ♦ ^ I D z 1 ] 

(32) 

For developing the explicit function of s, 
(A(s)] t corresponding to [A(ik>], the [D 0 ] and 
|D 2 ] contribution to [A(ik)J are identified 
separately, and the explicit occurrence of s in 
Equation (30) is maintained. 

[ A( ik,s) ] - |ZP] tNID(lk)]- 1 

* |D e ]. | |ZP] (NID(ik)]- 1 [D z) 


Let : 


lA^ik)] - IZP] [NID(ik)]- 1 [D 0 J 


(33) 

(34) 


Then : 


[ A( i k, s ) ] - [ A ? ( ik) | 4 | ( A z ( ik) ] 


(35) 


(36) 


Preliminary to approximating (A(ik,s)) by an 
explicit function of only s, [A x (ik)| and [A z (ik)] 
are approximated by [ A T ( p ) ] and [A z (p)] ( vhere p 
is the nondimensional form of s: p - cs/V. 

Following Reference [8}, the following terms 
are approximately, 


[A x (p)l - [B X q] 4 

n 

I 

j-1 

Pxjjp 

p 4 b . 
F 3 

(37) 

1 A 2 (p) 1 - lB z0 ] 4 

n 

t 

j-1 

Pzijp 

(38) 

P 4 bj 

These matrices can 

be 

obtained 

by generating 

aerodynamic matrices 

for 

several 

values of k and 

then employing a least-squares fit 



Because the state-space equation will be 
written in terms of s, Equations (37) and (38) are 
written in terms of s by letting p = cs/V: 


n [Bjj ] 

[A T (S>] - [B X0 ] 4 s l rTs 

j-1 

|A z (s)| - (B z0 l ♦ S I ftlL. 

j*l J 


(39) 


(40) 


vhere 

pj - v b j /<: 

Combination of Equations (31), (35), (36), 

(39), and (40) leads to the following approximate 
expressions for the aerodynamic forces: 


(Zaero) 


I* 2 


(B T o) 


P4 

j.i * * ‘i 


(z ) 


♦ ]■ pV 2 


, B . + I Pill 
i B zoi v s ♦ e 

j-1 


{*} 

(41) 
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5. INTEGRATED TOTAL BODY AND ELEMENTAL 


BODY DYNAMICS 

The general form of the linear flexible body 
equation is 


M'q + D q + G'(u>)q ♦ K'(oxI>,K)q 
+ A'(co,u>)q ^ L' ( u>, 5>, f s , p s ) 


(42) 


where q- [ u^, . , uj , , . . , uj , 0 ^ ]T M , D, and K 
are mass, damping, and stiffness matrices 
respectively of the airplane which are obtained by 
a traditional finite element method, such as 
NASTRAN. The other terms, G', K' , A' and L' are 
obtained from equations (16), and (42). Note that 
K' is a symmetric matrix. Equations (3), (11), 
(12), (13) and (42) can be written in the state 
space like format as follows: 


■M Tot e o o o' 


■V ■ 

o i* o r 0 


to 

0 0 E n 0 



MI EO e \| 0 M\ 


i • 

ir\ 2 \ 


o 

o 

o 

o 


V ' 


' F ' 

o ai* i r l 


0) 


P 

0 0 0 -E n 


*1 


0 

. 0 Ai K'+A' G' + D_ 


Ln 2 . 


. X . 


(43) 


r l ^2 3 (t”> s (2r s T u s E-r s u|-u s r s T ) 

* I ®s I s* II s®s) 01 ♦In s [(»r s ) - 

♦ r s «|u s . rai s p s 

E = 3x3 unit matrix 


( 48 ) 


H . block diagonal 6 nx 6 n matrix where block 
diagonals are 3x3 matrices 

- Block diagonal [ mjE, I lp . m s E, I s , 

. . , m n E, I n 1 


E e0 - l E 0 

E 0 . 

. E 0 1 T 

6nx3 matrix 

Eqe - [ 0 E 

0 E . 

. 0 E | T 

6nx3 matrix 

M' - M(E n - 

r E0 cT E0 M/M Tot) 

(49) 

Aq - M(E 0E - 

R) 


(50) 


r l 

0 

f 2 
0 

0 J 

A 1 = M C 0E R 


(51) 


(52) 


Equations (44)-(55) explain various terms In 
equation (43). 


Hi 

n 2 


(44) 

(45) 


0 is a null matrix and dimensions are context 
dependent . 


E n is 6 nx 6 n matrix. 

r 0 h 2 - lm s l s u 3 +tl s & 3 

E Tij - «(Em g (2r 3 T u s E-r s u|-u s r s T ) 

4 E0 s X s -n s 0 s ) (Em s (2r s T u s E 
- fs u I- u s r s T ) + 2: 0s I s- n s^s) w 


(46) 


(47) 


[ Eo£«p =» block diagonal matrix of dimension 
6 nx 6 n. Each block is 3x3. 

The diagonal blocks are 

[ w 0 « 0 . . u 0 ] 

I e EO w J~ ■ block diagonal matrix of dimension 
6nx6n. Each block is 3x3. 

The diagonal blocks are 

[0 & 0 « . . 0 u> ] 

lHEQgu>r - block diagonal matrix of dimension 
6nx6n. Each block is 3x3. 

The diagonal blocks are 

[ 0 0 (I 2 U)“ .. 0 (I n w)-] 

K' - K ♦ |Eoe w 1~ M [ Eq£«] ~ ♦ M { [EgowJ* 

* 1%)“!' -Ieo^eS M/M Tot ) (53) 
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A' - M[ £<)£«)” - IMEQE^r 

♦ H ( [Ieow]" - ^EO^EO^Tot ) 

♦ (EOE^T I mi: oe w 1'’ 


(54) 


C' - MI£o E w]" + [EoE w i~ M ~ l ME OE w l~ 

♦ 2H { [EeO w 1*' " E EO &E EO H/H Tot ) < 55 > 


Definitions of I Coqu] “ * [ ^EO 40 1 * an( ^ l M ^OE w l 
are very similar to [£qew]~, I ^£0^1 » an ^ ( M ^OE w i 
and hence they are not given here. 


Equation (43) has the form 


A 0 X + Aj X - U 


(56) 


where definitions 
Equation (56) can 



of Aq, , 
be written 


X 



X, and Q are obvious, 
as 


(57) 


Eq.(57) can be simply written as 
X = A X + B U 


(50) 


Definitions of A and B are obvious. 

Eq . (58) can be used for the time simulation. 
To better understand the interaction between the 
total body and the flexible body dynamics, steady 
state maneuvers (i.e. constant angular rates of 
the body reference frame) are studied. 


By putting the derivative of w to zero and 
including the aerodynamic force representation 
from equation (41), equation (43) becomes: 


I 

-pV B 
0 
0 


~r 


zo 


0 

-pV B 
I 
0 


~T~ 


Z1 


0 

0 

I 




0 -I 

K'*A'-pV 2 B_ n G' + D' 

~y T0 

0 -I 

0 -I 



Q 

0 

0 

0 


( 59 ) 


Note that the coefficient matrices on LBS of 
the equation (59) are time invariant when the 
angular rate, w, is constant. Hence the 
eigenvalues of the system can be used to check the 
stability of the system and to study the effects 
on modes of the system at different angular rates. 

6. APPLICATION TO FI KITE ELEMENT HODEL 

A large-order finite element model (FEM) of a 
generic fighter was obtained for use in the 
application of this method. The aircraft planform 
is similar to an F/A-10, although stiffness and 
mass data do not necessarily represent this 
airplane. Although the FEM consists primarily of 
beam elements, it is a highly detailed model 
containing an A-set of 228 degrees of freedom 
(DOF) and approximately 200 structural elements. 
Aerodynamic modeling of the aircraft consisted of 
230 boxes, and can be seen in Figure 2. The 
doublet lattice method was used to formulate 
aerodynamic influence coefficient matrices. Eight 
values of reduced frequency were used to calculate 
unsteady aerodynamic matrices. 

Certain assumptions used to develop the 
equations required that some modifications be made 
to the model. The equations assume that the 
mathematical model has six DOF for every subbody. 
If these matrices are generated from a FEM, this 
is rarely true. In NASTRAN , this corresponds to 
the initial global set (G-set) of coordinates. 
These DOF cannot normally be used, however, 
because many are constrained due to the method of 
modeling and imposition of boundary conditions. 



Figure 2. Aerodynamic Configuration 
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Table 3. Flexible Mode Frequency and Damping for Roll Rate Maneuvers 
Full Stiffness 


msc/n; 

FLUTl 

ANAL' 

^STRAN 

rER 

fSIS 

FUG 
ANAL 
0.0 D 
ROLL 

MR 
fSIS 
IQ / SEC 
RATE 

FLXAIR 
ANALYSIS 
90.0 DEC/SEC 
ROLL RATE 

FLXAIR 
ANALYSIS 
180.0 DEG/SEC 
ROLL RATE 

FLXAIR 
ANALYSIS 
240.0 DEG/SEC 
ROLL RATE 

FREQUENCY 

Hz 

DAMPING 

FREQUENCY 

Hz 

DAMPING 

FREQUENCY 

Hz 

DAMPING 

FREQUENCY 

Hz 

DAMPING 

FREQUENCY 

Hz 

DANPINC 

6.734 

B.948 

9.085 

14.192 

16.434 

18.736 

21.172 

23.172 
24.352 
29.578 
32.916 

-.0689 

-.000459 

-.0740 

-0.175 

-.0730 

-.0249 

-.0188 

-.0268 

-.0406 

-.00707 

-.0458 

6.699 
8.957 
9.056 
14.371 
16.779 
18.756 
21.812 
1 23.333 
24.692 
29.719 
33.439 

-.0653 

-.000422 

-.0722 

-0.142 

-.0806 

-.0297 

-.0157 

-.0292 

-.0416 

-.00878 

-.0387 

6.696 

8.955 

9.053 

14.368 

16.778 

18.755 

21.813 

23.333 

24.691 

29.276 

33.448 

-.0653 

-.000467 

-.0722 

-0.142 

-.0806 

-.0298 

-.0159 

-.0293 

-.0416 

-.00871 

-.0387 

6.685 

8.946 

9.044 

14.361 

16.774 

18.753 

21.814 

23.331 

24.689 

29.721 

33.437 

-.0655 

-.000539 

-.0722 

-0.142 

-.0806 

-.0298 

-.0165 

-.0294 

-.0416 

-.00865 

-.0387 

6.673 

8.936 

9.035 

14.353 

16.771 

18.751 

21.815 

23.333 

24.688 

29.722 

33.436 

-.0656 
-.000602 
-.0722 
-0. 142 
-.0805 
-.0298 
-.0162 
-.0294 
-.0416 
-.00863 
-.0387 


These constrained DOF present a problem which 
requires either the modification of the equations 
or of the input matrices. 

Another assumption made in the equations is 
that the mass matrix is block diagonal. However, 
the typical mass matrix from a FEM analysis 
contains coupling terms. These arise because of 
the following reasons: 

1) Mass data may be input at locations other 
than structural grid point locations. 

2) Coupling results from the use of 

dependency relations (multi-point 

constraints in NASTRAN). 

3) Coupling results from the static reduction 
if inertia is lumped on any of the 
omitted DOF (Guyan reduction). 

These considerations make it necessary to 
adjust the model as follows: 

1) The inertia is relumped so that it is 
located at exact grid point locations. 

2) Inertia at dependent DOF is relumped so 

that it is associated only with 

independent DOF. 

3) Inertia located at DOF which are 

eliminated by the Guyan reduction process 
must be relumped at retained DOF (A-set). 


4) A Boolean transformation matrix is formed 
for use in expansion of the FEM A-set DOF 
to the 6n DOF required by the equations. 
After forming the state space equations, 
this same matrix can be used to eliminate 
those DOF. 

A NASTRAN flutter analysis of the vehicle was 
conducted for a case representing Mach .7 and an 
altitude of 20,000 ft. Mass, stiffness, and 
aerodynamic matrix data were obtained from NASTRAN 
for this case. The necessary matrices for the 
rational function approximation of the 
aerodynamics were obtained by a least squares fit 
using aerodynamic matrices for reduced frequencies 
of 0.0, 0.2, and 0.8 . The state space equations 

were formed and eigenvalue solutions were obtained 
for various values of roll rate. 

For zero roll rate, the results agreed with 
the NASTRAN analysis. Increasing roll rate shoved 
little effect on the flexible modes of the system, 
as can be seen in Table 3. The rigid body modes 
were affected, however. A root locus plot of the 
rigid body roots as a function of roll rate is 
shown in Figure 3. For zero roll rate, two stable 
real roots and one stable complex conjugate pair 
are obtained - corresponding to a roll convergence 
mode, a spiral mode, and an oscillatory dutch roll 
mode. With increasing roll rate, however, we see 
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Figure 3. Rigid Body Eigenvalues for Roll Rate 
Maneuvers 
Full Stiffness 



REAL AXIS 


LEGEND 

o o deg/sec 

A 10 DEC/SEC 
□ 20 OEC/SEC 

X 30 DEC/SEC 
o 40 DEC/SEC 
^ 60 OEC/SEC 

H 70 DEC/SEC 
■ 80 0 EG/SEC 

♦ 90 DEC/SEC 

© 100 DEC/SEC 

3 HO 0 EG/SEC 
ffl 120 DEG/SEC 
® 150 OEG/SEC 

E3 180 DEG/SEC 

• 210 OEG/SEC 

O 240 DEG/SEC 


that some roots become unstable, and also change 
from real to complex and back again to real. 
Another case, representing a more flexible 
airplane, shows the same behavior (Figure 4), 
although the changes occur at lower roll rates. 
This case represents 50 X of the initial overall 
airplane stiffness. Table 4 shows again that the 
flexible modes were not greatly affected, even for 
the reduced stiffness case. 

7. DISCUSSION 

Dynamic equations have been derived for a 
flexible fixed wing aircraft, including an 
explicit representation of unsteady aerodynamic 
forces. The aircraft is assumed to be a 
collection of elastically interconnected discrete 
rigid subbodies. Deformations are assumed to be 
small in the body reference frame, thus allowing 
the equations to be linearized in the deformation 
variables. 


The hybrid coordinate method is used to 
derive the total body and the elemental body 
dynamic equations which are then converted to 
matrix form. These equations are integrated in a 
state space format, along with a rational function 
approximation of the unsteady aerodynamic forces. 
These equations can be used for simulation. For 
the case of constant angular velocities of the 
body reference frame, the coefficient matrices 
become time invariant, allowing the use of an 
eigenvalue analysis to evaluate the effects of the 
angular rates on the system dynamic properties. 
Vhen this method is applied to a realistic finite 
element model of a generic high-performance 
fighter, significant changes in the stability 
characteristics of the aircraft are observed. 
With increasing roll rate, some roots become 
unstable, and also change back and forth from 
complex to real. The dutch roll mode becomes two 
real roots, one of which combines with the spiral 
mode to produce an unstable oscillatory mode. The 
other real root from the original dutch roll mode 
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Table 4. Flexible Mode Frequency and Damping for Roll Rate Maneuvers 
50 Percent Stiffness 


FLXAIR 
ANALYSIS 
0.0 DEG/SEC 
ROLL RATE 

FLXAIR 
ANALYSIS 
60.0 DEG/SEC 
ROLL RATE 

FLXAIR 
ANALYSIS 
90.0 DEG/SEC 
ROLL RATE 

FLXAIR 
ANALYSIS 
180.0 DEG/SEC 
ROLL RATE 

FREQUENCY 

DAMPING 

FREQUENCY 

DAMPING 

FREQUENCY 

DAMPING 

FREQUENCY 

DAMPING 

Hz 


Hz 


Hz 


Hz 


4.860 

-.07367 

4.858 

-.0737 

4.855 

-.0737 

4.839 

-.0739 

6.335 

-.000492 

6.333 

-.000553 

6.331 

-.000590 

6.318 

-.000727 

6.500 

-.0929 

6.500 

-.0930 

6.500 

-.0930 

6.484 

-.0930 

10.230 

-0.186 

10.229 

-0.186 

10.227 

-0.186 

10.216 

-0.185 

11.871 

-.0073 

11.871 

-.0874 

11.870 

-.0874 

11.867 

-.0873 

13.334 

-.0431 

13.334 

-.0431 

13.333 

-.0432 

13.331 

-.0431 

15.435 

-.0196 

15.435 

-.0199 

15.435 

-.0200 

15.436 

-.0204 

16.527 

-.0364 

16.527 

-.0365 

16.527 

-.0366 

16.525 

-.0367 

17.498 

-.0578 

17.497 

-.0578 

17.497 

-.0578 

17.494 

-.0578 

21.021 

-.0125 

21.021 

-.0124 

21.022 

-.0124 

21.024 

-.0122 

23.671 

-.0479 

23.711 

: -.0479 

23.671 

-.0479 

23.669 

-.0479 


Figure 4. Rigid Body Eigenvalues for Roll Rate 
Maneuvers 
50 Percent Stiffness 
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ns can be observed for the 50X stiffness case. 


combines with the roll convergence to form another 
oscillatory mode vhich becomes more stable with 
increasing roll rate. The effect on the flexible 
modes of the aircraft was minimal for this 
configuration. The behavior of the rigid body 
modes is somewhat dependent on airframe stiffness, 


It is expected that a design with increased 
span would show a greater effect due to roll rate 
for both the rigid body and flexible modes. This 
should be given consideration in the design of any 
future high-performance aircraft. 
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